Tunneling of Bloch electrons through vacuum barrier 
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Tunneling of Bloch electrons through a vacuum barrier introduces new physical effects in com- 
parison with the textbook case of free (plane wave) electrons. For the latter, the exponential decay 
rate in the vacuum is minimal for electrons with the parallel component of momentum ky = 0, and 
the prefactor is defined by the electron momentum component in the normal to the surface direction. 
However, the decay rate of Bloch electrons may be minimal at an arbitrary ky ("hot spots"), and 
the prefactor is determined by the electron's group velocity, rather than by its quasimomentum. 



There is a general feeling in the applied physics com- 
munity nowadays that the next decade will bespeak 
an advent of magnetoelectronics n exploiting spin, rather 
than charge, degrees of freedorcu. Many of such spin- 
tronic devices are based on the phenomenon of quantum 
tunneling, and specifically on the difference between tun- 
neling currents in different spin channels, this opens up a 
possibility to control the electric properties via magnetic 
field. In view of all that, an explosion of publication on 
spin-polarized tunneling which started around 1995, is 
not surprising at all. 

Interestingly, despite the fact that the tunnelling is one 
of the best studied phenomena in quantum mechanics, 
there is still substantial diversity in microscopic under- 
standing of tunneling in real systems. Tunneling prob- 
lems are very easily solved in one dimension and for free 
electrons; but it is not so obvious how this is related to 
real systems, and how to incorporate the effects of real- 
istic electronic structure. This is the reason why most 
theoretical papers use the free-electron model, that is, 
the electronic wave functions that are plane waves and a 
spherical Fermi surface. The deviations of the wave func- 
tion from the single plane wave form, and of the Fermi 
surface from a sphere, are very often crucial for under- 
standing the physics of tunneling. In particular, a totally 
counterintuitive result has been observed in some recent 
calculations!], when the electrons with nonzero quasimo- 



mentum parallel to the interface had a larger probablility 
to tunnel through a vacuum barrier compared to those 
with zero parallel quasimomentum. It was also pointed 
outtj that electrons in different Bloch states with the same 
energy and quasimomentum may, in principle, have dif- 
ferent decay rates in vacuum: another counterintuitive 
result. In this regards, it is important to establish a 
formal theory for tunneling of Bloch electrons through a 
vacuum barrier, elucidating the qualitatively new aspects 
of this process as opposed to the free electron tunneling, 
particularly because only the latter is discussed in the 
classical textbooks. This is the goal of the current paper. 

The simplest ease of a tunneling contact is the so-calles 
Sharvin contacts, which is essentially an orifice between 
two metals (or a metal and the vacuum), whose size is 
smaller than the mean free path of electrons in the bulk. 
All electrons with the a positive projection onto the cur- 
rent direction (which we will denote as x) pass through 
the contact. Conductance of a Sharvin contact between 
two identical metals is 
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where A is the contact area, N is the volume density of 
electronic states at the Fermi level, v is the Fermi veloc- 
ity, and brackets denote Fermi surface averaging: 
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Integration and summations are over tk&,states with 
v\tia,x > 0, and fl is the unit cell volumetru. k, i, and 
a denote the quasimomentum, the band index, and the 
spin of an electron, respectively. This formula can be 
derived by considering the voltage-induced shift of the 
Fermi surfaceQ, but there is a more instructive derivation 
starting from the Landauer-Buttiker formula for the con- 
ductance of a single ballistic electron, Go = e 2 /h. In this 
formalism, the total conductance is equal to Go times 
the number of conductivity channels, N CCl which is de- 
fined as the number of electrons that can pass through 
the contact. Assuming that the translational symmetry 
in the interface plane is not violated, we observe that 
the quasimomentum in this plane, ky , is conserved, and 



N cc is the number of quantum-mechanically allowed ky 's. 
Thus N cc is given by the total area of the contact times 
the density of the two-dimensional quasimomenta. The 
latter is simply S x / (2tt) 2 , where S x is the area of the pro- 
jection of the bulk Fermi surface onto the contact plane. 
Thus 
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This is an important result. To the best of our knowledge, 
Walter Harrison was the first to spell it out in 19610, 
and there is no lack of more recent papers manifesting 
proper understanding of this issue (e.g., Ref.El). How- 
ever, till now many otherwise correct and useful papers 
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erroneously identify the number of conductivity channels 

incorrect] 

Eq. H is the basis for all more sophisticated expres- 
sions describing various aspects of quantum tunneling. 
None of them may explicitely depend on the bulk den- 
sity of states. It may, however, be that that instead of the 
straight (A^i^l) averaging one has to compute a weighted 
average, with the weights coming from tunneling matrix 
elements, or other additional physics. 

Eq. [j] takes care of one important difference between 
the free electrons and the Bloch electrons: deviation of 
the Fermi surface from a sphere, for S x 7^ irkp. Another 
important difference that is often neglected is that be- 



and the density of states at the Fermi level, that is, 




tween the group velocity K~ de^/dh and the phase veloc- 
ity fik/rriQ. One can get some qualitative understanding 
of the role that this fact plays in tunneling by considering 
a simplified model, where electrons in metal are approx- 
imated by free electrons with an effective mass different 
from the free electron mass. In this approximation, the 
effective mass is responsible for the difference between 
(Kk./mo) and (de^/hdk), and the transparency of a sym- 
metric rectangular barrier is defined by the standard for- 
mula (see, e.g., Ref.B): 
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where wl(_r) stands for (k-dependent) Fermi velocity in 
the left and in the right leads, and the imaginary quasi- 
momentum fiK is calculated from the energy conserva- 
tion condition, 



U + h 2 [k\\ 2 



K 2 ]/2m = E, 



where mo is the free electron mass, U is the barrier 
height, and d is its thickness. 

The physical reason that one has to use group veloc- 
ities, and not wave vectors, is very profound and ex- 
tends well beyound the limited scope of the effective 
mass model: these factors appear in the EqJ^ as a re- 
sult of matching the gradients of the wave functions at 
the interface, and the gradient is, in fact, the velocity 
operator for the Bloch waves. Another way to express 
the same idea is to recall the physical meaning of the 
usual quantum-mechanical requirement that the wave 
functions be smooth: it is needed to ensure the flux conti- 
nuity and, therefore, particle conservation. On the other 
hand, the expression for K includes the momentum, kii„ 
and the free electron mass, mo, because it comes from the 
solution p£ the Schrodinger equation inside the barrier (in 
vacuum)H3. The conductance of a contact described by 
Eq. |5| is given by the appropriately modified Eq.|I|: 

G= TSE^- S ^ ], W. (6) 

k 

It is instructive to consider the last formula in some 
limiting cases. First, let us consider a specular barrier. 
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It is defined by the limit U 
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co, d -> 0, Ud = V. Then 
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Note that in the literature the ratio V/hv x — Z is com- 
monly used to characterize the barrier strength. In prin- 
ciple, this quantity is different for different electrons, as 
v x depends on k. In the limit of low transparency, Z ^> I , 
D(k) = h 2 v F v R /V 2 . Substituting this into Eq.pl, we find 
that the total current is proportional to 
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where summation is, of course, over those k that are 
allowed in both left and right lead. Roughly speak- 
ing, the total conductance is defined by the smaller of 
the two (NvD's, that is, by mm ({Nvl) L , (Nv 2 ) R ) . In 
the high transparency limit D is still smaller than 1, 
D = 4vlv r /(vl + v R ) 2 , (so-called Fermi velocity mis- 
match), but in most cases this is not a large effect: factor 
of two mismatch reduces D by only 10%. 

In the case of a thick barrier, defined as dK 3> 1, Eq. 
H can be expanded in h 2 k 2 /4mo(U — Ep), and the trans- 
parency is 



{U-E F + m v 2 /2){U - E F 
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where y/2rrio(U - 
proportional to 



Eo)/hd = W <C k (thick barrier limit). W does not depend on k. The tunneling current is 
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Let k n be the set of points on the Fermi surface where k|| = (note that for Bloch electrons beyond the effective 
mass approximation tunneling from some of these points may be suppressed by symmetry, as discussed later in the 
paper). Except in the exponent, we can put k || to zero, 

T 1 \- \ fjiu r_*!l\ 2m 2 {U - E F )v L v R 

J0C (2^2^17 fc6Xp[ W l j (U-E F + m vi/2)(U-E F + m v%/2) [n) 

VL VR 



E 



m h 2 vl/2 + U - E F m h 2 v R /2 + U - E F 



All omitted in this expression factors are k ||- 
independent. One should not be confused by the fact 
that, unlike Eq. ||, the numerator here does not have the 
third velocity. We have reduced our problem to an effec- 
tive ID problem, in which case the role of the density of 
states is played by the inverse velocity. Correspondingly, 
the product Nv cancels out. 

Eqs. ||,^] emphasize the role of kinematics in tunneling. 
For instance, the long-standing problem of the reversed 
(compared to the density of states) spin polarization of 
the 3d ferromagnets is entirely explained in terms of kine- 
matics. Direct calculations show that s-like electrons in 
Fe, Co and Ni have much larger Fermi velocity than <i-like 
electrons. Taking this fact into account brings the cal- 
culated spin polarization to a very good agreement with 
experiment, without making any additional assumptions 
about the character of the surface statesBO. This is by no 
mean surprising: the bulk transport is controlled by the 
same factor (Nv' 2 .') , and the Ohmic current in these met- 
als is carried predominantly by s-like electrons. It is only 
natural that in another transport phenomenon, tunnel- 
ing, these electrons also play the leading role. We would 
like to emphasize that the effect considered above (as op- 
posed to another effect discussed later in the paper) is not 
related to the s or d symmetry of the wave functions, but 
to the group velocities in the respective bands. In other 
cases the "light" and the "heavy" bands may not be di- 
rectly related to the angular symmetry of the wave func- 
tions. For example, in SrRu03 both spin-up and spin- 
down Fermi surfaces are made up by Ru t2 g d-electrons, 
but the average group velocity in the spin-majority cha»- 
nel is twice smaller than that in the spin-minority one£3. 
As a result, although the spin polarization of the density 
of states is positive, N-\ > while the_transport spin 
polarization is negative, {Nv}-\ < (Nv)±E3 

Now we have some understanding of the two remark- 
able differences between the free electrons and the Bloch 
electrons: the effect of the Fermi surface geometry and 
the difference between the group and the phase velocities. 
There is, however, yet another, extremely important, dis- 
similarity between the two systems, recently pointed out 
by W. Butlem: the difference between the momentum 
and the quasimomcntum. In order to discuss this dif- 



ference, and its physical consequences, let us consider 
reflection of an individual Bloch wave from a metal sur- 
face. Let x be the direction normal to the surface, and 
r the coordinate in the surface plane. At x < we have 
a metal, and vacuum at x > 0. The vacuum potential is 
again U, and the Fermi energy is E. Since we have perfect 
in-plane periodicity, the wave function at any x can be 
classified by k||, and is given by 

^(k||,a;,r||)=^exp[i(k||+G)r||]F G (k||,a:). (12) 

G 

The quasimomentum in the surface plane, Tjkii , is con- 
served, as well as the energy. In vacuum, the solution of 
the Schrodinger equation is 

^ T (k||,a:,r||) = ^a G exp[i(k|| + G)r||] exp(-K G x), 

G 

(13) 

where G is the 2D reciprocal lattice vector, and Kg is 
now defined taking into account the kinetic energy associ- 
ated with the given reciprocal lattice vector, U + h 2 [h\ 2 ~ 
K G ]/2m = E. An incoming Bloch wave with a given k 
penetrates into the barrier as a linear combination ( |l3| ) 
with the coefficients oq defined by matching conditions, 
set by the requirement of continuity of the wave function 
and its derivative: 

f g(0) exp[i(ky + G)r,|] =^a G exp[t(k|| + G)r,|] 

G G 

*g(0) exp[i(k|| + G)r,|] = - £ a G K G exp[*(k u + G)r,|]. 

G G 

since this has to hold for any rn, a G = -Fg(O) , and 
F' G (Qi) = -a G K G for each G. Thus 

F G (0)A- G + F G (0) = 0. (14) 

If F G (x) were a linear combination of the bulk Bloch 
waves with the energy E and the quasimomentum in the 
plane Tikii , 

F G (x) = Uk x (x) exp(ik x x) + au_ fex (a;) exp(-ik x x), (15) 



3 



we would have only one free parameter, a, to satisfy Eq.|l4| 
for all G's, which is obviously impossible. The answer is 
that Fg (x) has the form ( |l5| ) only far away from the sur- 
face, while near the surface it is distorted as required by 
Eq.|l4]. This emphasizes once again the role of of surface 
states in tunneling. In fact, one of the ways to realize 
the necessity of forming the surface states is that the 
bulk Bloch functions, in general, cannot be augmented 
continuously and smoothly into vacuum. 

In the case of a thick barrier, the actual tunneling cur- 
rent will be defined by that component of the wave func- 
tion ( |l3| ) which has the smallest K, that is, by the one 
with G = 0. The amplitude of this evanescent wave is set 
by ao- As pointed out by Butlcrcl, k|| = is a high sym- 
metry direction (TX), and the electronic states possess 
certain symmetry in the yz plane. In particular, ao for 
some states may vanish by symmetry, in which case the 
decay rate K will be defined by the smallest G allowed 
by symmetry. Since we consider now a thick barrier, this 
essentially means that tunnelling from such a band will 
be defined not by the k|| = state, but, rather counter- 
intuitively, by general (not high symmetry) points in tic 
2D Brillouin zone (as confirmed by actual calculationsa) . 
Indeed, consider a band where by symmetry Fq(0,x) = 
at k|| =0. At k|| ^0 thus F (k\ l: x) = F£(0,x)k 2 , 

while K = ^2m (U - E)/h 2 + k 2 » ^/2m {U ~ E)/h+ 

hkj/2- s /2m Q (U - E) = K + kj/2K . The optimal dis- 
tance from the zone center that gives maximal contribu- 
tion to the tunneling current can be estimated by maxi- 
mizing with respect to fc|| of 

F^'{0,x)k 2 eiq>(-K d - kfd/2K ), 

where d is the barrier thickness, which gives fc|| ~ 
y/2Ko/d. For Fe, for instance, Kq « 0.6 a.u., about the 
same as the TX distance. Thus for a barrier, say, of 5 lat- 
tice parameters, fc|| ~ 0.2 a.u., a sizeable distance from 
the center of the Brillouin zone. Yet another counter- 
intuitive result is that the low transparency limit is not 
unique: in the thick barrier limit tunneling is predomi- 
nantly from the states infinitely close to the zone center. 
However, in the high barrier limit, which is another way 
to implement a low transparency asymptotics, tunneling 
occurs far away from the zone center, possibly at the zone 
boundary. This is the effect observed in Refs.13. 

To conclude, we discussed here three new effects which 
appear in tunneling of the Bloch electrons through a vac- 
uum barrier, as compared with the textbook case of free 
electron (plane wave) tunneling. These effects are due 



to (i) complexity of the Fermi surface geometry ( "fermi- 
ology"), (ii) difference between the group and the phase 
velocities of a Bloch electron, and (iii) nonconservation of 
the parallel component of electron momentum (and con- 
servation of its quasimomcntum) . Each effect influences 
the tunneling current in its own way, and as a result even 
for the most simple case of a vacuum barrier, the tun- 
neling of the Bloch electrons appears to be qualitatively 
different from the free electron tunneling. 
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